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' This manuscript reports a stochastic dynamical scenario whose associated stationary probability 

density function is exactly a previously proposed one to adjust high-frequency traded volume dis- 
tributions. This dynamical conjecture, physically connected to superstatiscs, which is intimately 
related with the current nonextensive statistical mechanics framework, is based on the idea of local 
ft , fluctuations in the mean traded volume associated to financial markets agents herding behaviour. 

q_) . The corroboration of this mesoscopic model is done by modelising NASDAQ 1 and 2 minute stock 

Pj^ ' market traded volume. 
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^ \ In the last years we have assisted to a rapidly growing interesting of physicists in the analysis of financial systems 
g ■ and the description of their statistical and dynamical properties Actually, some of the so-called stylised 
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facts such as non-Gaussianity in relative price change (the return) probability density functions (PDFs) 0, 0,M El 
flfl multi- fractal behaviour of return time series ll 2L ll 3l Hij. asymptotic power-law like auto-correlation funtion 
in the standard deviation (the volatility) of return [15l llfij. among many others, were (also) verified in a physical 
context. Besides their own relevance, these empirical observations were important to the establishment of either 
new theoretical models, which mimic reality much better than previous ones [la. flTL llSl Il9|. or connections between 
financial systems and already existing statistical frameworks [ToL I20L [2lj ]. particularly, the current non-extensive 
1 I statistical mechanics |p] , based on Tsallis entropy |2^ . which, e.g., allowed the derivation of a new and closed model 
for option-pricing [ HTmI^ . 

Additionally to macroscopic observables as the return and the volatility, mentioned above, another significant 
quantity in financial markets is the traded volume, V, which corresponds to the number of shares traded (in a 
determined period of time, T). This is, in fact, along with volatility, an important parameter in trading strategies 
assumed by many brokers |1EI23- I n fig-fflis presented, as mere illustration, the daily normalised traded volume of 
Microsoft Corp. in NASDAQ exchange from 2000 to the beginning of 2005. 

In a previous work, P. Gopikrishnan et al. j28j reported the evidence, for large volumes, of a power-law behaviour 
fT) I in the stationary PDF. The analysis of traded volume stationary PDF was later extended to all range of values by R. 
Osorio, L. Borland and C. Tsallis (OBT) who proposed a fitting expression that presents an amazing precision. 
The respective formula is, 

1 f v\ a ( v\ 

" l— ! ! 

where v represents the traded volume expressed in its mean value unit ((V)), i.e., v t = -^y Z — J °° (j) a exp ? (— f ) dv, 
. at and 9 are positive parameters and 

'2: ^ 

exp (x) = [1 + (1 - q) x] i-« , (ex Pl (x) = e x ), (2) 

o . 

^ ' is the q-exponential function, which emerged in the context of the non-extensive statistical mechanics 0. 

In the sequel of this article, through a set of liable arguments, we will introduce a stochastic dynamical mechanism, 
which has eq. JQ) as stationary solution for its probability density function and then compare both analytical and 
simulation results using values obtained from NASDAQ stock market. 



Consider the following stochastic differential equation, 



dr = - I v - 2L+_J dt+J2vl dW t , (3) 
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FIG. 1: Evolution of Microsoft Corporation daily normalised traded volume in NASDAQ exchange from the 21 s January 2001 
to the 21 s * January 2005. Each tick represents 6 months. 



(Wt represents a regular Wiener process following a normal distribution with null mean and unitary variance) as the 
traded volume updating law 1 . As can be seen, eq. © is composed by a deterministic term, similar to a restoring force 
and a stochastic term, which is a typical example of multiplicative noise very usual in problems presenting clustering 
profiles, like the one treated in this manuscript (see fig. Q), intermittence j^ij or turbulence |3ll Is^. The first term 
can be interpreted as a restoring market mechanism (with constant 7), which aims to keep the traded volume of the 
stock at some typical value = ^-jjp- . 2 The second one reflects stochastic memory, and mainly, the upshot of large 
traded volumes. In effect, large traded volumes will induce large values in the stochastic term, but, due to Wt random 
character(sign), it could lead to an increase (W t > 0) or a decrease (W t < 0) of v, i.e., to a stirred or serene state of 
the stock. 

For eq. J2J) we can write the corresponding Fokker-Planck equation [3j| , 



dp(v, t) d 
dt dx 
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p(v,t) 



dx 2 



v-p(v,t) 



(4) 



and determining its stationary solution we obtain 



p(v) = 



r[i 



(Pvf exp(-0v) (a>-l, /3>0), 



(5) 



which is just the Gamma distribution, with mean value, (v) = ^-jp, standard deviation, (v — (v)) 2 = and most 



probable value, v* = % ( 



dv 



0). 



Assume now that, instead of constant in time, is now randomly varying on a time scale much larger than the time 
scale of order 7 -1 necessary by eq. Q to reach stationarity, hypothesis that is equivalent to consider local fluctuations 
in the mean value of v |34j , akin to the proposal of moving average in the volatility |35l | . These fluctuations can be 
microscopically associated to changes in the volume of activity defined as the number of traders who performed 
transactions of the stock at that time [3(|. Furthermore, let us assume that f3 follows a Gamma (stationary) PDF, 



1 



XT [5] 



s-i 



exp 



(5 > 0, A > 0). 



(6) 



This assumption transforms eq. (0) meaning from the absolute probability of v into the conditional probability of v 
given P, 



p(v) ->p(v[ 



(3 



r[i 



(flv) a exp (—(3v) 



(7) 



1 Eq. i s i n f ac t defined in mathematics as Feller process l2El . According to its definition the traded volume V is an integer quantity. 
However, the normalised traded volume is a rational quantity and so perfectly compatible with a continuous approach. 

2 Following v definition, should equals 1. Nonetheless, we will keep it in the presented form for reasons that will become clear later on. 



It is noteworthy that the problem of macroscopic non-equilibirium systems with complex dynamics in stationary 
states with fluctuations of intensive quantities on long time scales and its connection with non-extensive statistical 
mechanics, based on entropy |22^ . 
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(Si = — y^.- pi In pi = Sbg where BG stands for Boltzmann-Gibbs) was firstly presented by G. Wilk and Z. 
Wlodarczyk [23], later applied by C. Beck to the case of turbulent fluids j^t and finally extended and defined 
by C. Beck and E.G.D. Cohen as super statistics pKl Hcj. 

For this case, we have that the joint probability of obtaining a certain value of v and certain value of /3, is 
P (v, (3) = p (v\P) P (/3) and the marginal probability of having some value v is given by 



P(v) 



P(v,j3) dp 



p(v\0) P(I3) dp. 



Using eq. 10 and eq. J7J in eq. (JSJ) and performing the integration we get, 

\r[l + a + 5] 



P(v) = 



Doing the change in parameters 



T[l + a] T[S] 



{Xv) a (1 + Xv) 



-l-a-S 



(a + 5> -I) 
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a — 1, 



(8) 



(9) 



(10) 



one can write eq. © as 



P(v) 



(q-l) a+1 T 
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er 




T[a + l] 



1 + (q - 1) ■ 



1/(1-9) 



or, using the g-exponential function definition, eq. 



P(v) 



(q-l) a+1 T 


1 




6T 


-±-r-a-l 
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(ii) 



which is exactly the expression introduced by OBT. Due to its structure, PDF (|lll) might be called q-generalised 
Gamma probability density function. The traditional Gamma PDF, eq. J2J, can be recovered from eq. (|11J) by 
considering q = 1 limit, which is equivalent to assume the inexistence of fluctuations in P and so P (P) equals the 
Dirac delta function centred in 6 . For small values of v, P (v) goes to zero as a power law, P (v) ~ v a , and for 
large values of v, it behaves asymptotically as another power law, P (v) ~ ti Q /(i-9)_ As a matter of fact, this is an 
appealing form, since it congregates, in a single expression, the two power laws used in financial practice. 

In order to verify the correctness of the result, we performed a set of numerical realisations, using an algorithm 
based on cq. and on a random generator for /3 whose associated probability is eq. ©. The values of parameters 
a, A and 5, presented in table [fl were obtained after fitting the 1 minute and 2 minute traded volume probability 
density functions of the ten high- volume stocks in the NASDAQ exchange during 2001, fig. El with eq. JJJ (equivalent 
to eq. dJ). The values of q and 9 are obtained using the relations presented in eq. IjlOl) . 

The numerical realisations for 1 minute and 2 minute traded volume are depicted in fig. [21 and fig. 01 
In fig. ^ fig. [SI and fig. 0] (Centre plots) it is clear the clustering behaviour, very similar to volatility time series, 
typical of a multiplicative noise dynamics. In fig. [21 and fig. 0] (Right plots), we can compare numerically obtained 
relative frequency, H (v), with analitycal results corresponding to values presented in table[U As it is straightforwardly 
visible, numerical and analytical results are in perfect accordance. Besides the possibility of mimic traded volume 
time series, this dynamical conjecture can certainly be used in other problems with time varying positive quantities 
associated to a PDF (QJ , such as granular material systems • 



To summarise, in this manuscript we presented a dynamical conjecture, based on a multiplicative noise mechanism 
and on the fluctuating character of some parameter of that system, which is able to provide the stationary probability 




FIG. 2: 1 minute (Left plot) and 2 minute (Right plot) traded volume probability density function P(v), for the ten high- volume 
stocks in NASDAQ exchange during 2001 vs. v. The symbols represent real data and lines p(v) from eq.Q, with q, a, values 
shown in table (After OBT in ref.[10]). The x 2 (per point) error function values are 1.8 x 10~ 3 and 6.5 x 10 -4 , respectively. 
For very low volumes is visible a small discrepancy, due to lack of small transactions (statistical effect). 

TABLE I: Fitting parameters obtained in ng|5]and dynamical parameters obtained through ea llOl 





q 


a 


9 


A 


<5 


1 min 


1.19 


0.93 


0.23 


0.826 


3.33 


2 min 


1.16 


1.36 


0.2 


0.8 


3.89 



density function previously introduced to shape the high-frequency traded volume stationary probability distribution 
of shares in stock markets. The present scenario, related with superstatistics and thus with non-extensive statistical 
mechanics, contains two ingredients that are often considered as elemental in the dynamics of such complex systems: 
memory (via multiplicative noise) and microscopic mechanisms of herding between traders (through the local temporal 
fluctuations of f3 or the v mean value) caused, e.g., by news, rumours or simply price movements. The description of a 
completely related return-volatility-volume mesoscopic dynamical theory, within non-extensive statistical mechanics 
framework, for financial markets is currently in progress. 
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FIG. 3: Numerical realisation for 1 minute traded volume of the ten high-volume stocks in NASDAQ exchange during 2001 
using the values of table [I] Left Plot: temporal evolution of (3 parameter; Centre Plot: temporal evolution of traded volume, 
v. Right Plot: Relative frequency H (v) vs. v. The symbols were obtained numerically and the line from JJJ^J P (v) dv with 
Av — 5 x 10~ 3 . The ratio between f3 updating time and stationarity time scale of v is 10 2 . 




FIG. 4: Numerical realisation for 2 minute traded volume of the ten high-volume stocks in NASDAQ exchange during 2001 
using the values of table [I] Left Plot: temporal evolution of (5 parameter; Centre Plot: temporal evolution of traded volume, 
v. Right Plot: Relative frequency H (v) vs. v. The symbols were obtained numerically and the line from JJJ^J P (v) dv with 
Av = 5x 10~ 3 . The ratio between [3 updating time and stationarity time scale of v is 10 2 . 
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